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Abstract. - We analyze the fidelity decay for a system of interacting bosons described by a Bose- 
Hubbard Hamiltonian. We find echoes associated with "non-universal" structures that dominate 
the energy landscape of the perturbation operator. Despite their classical origin, these echoes 
persist deep into the quantum (perturbative) regime and can be described by an improved random 
matrix modeling. In the opposite limit of strong perturbations (and high enough energies), classical 
considerations reveal the importance of self-trapping phenomena in the echo efficiency. 



While the physics of the last century was mainly charac- 
terized by great advances in understanding the properties 
of one-particle systems, recent experimental developments 
have put the effects of interacting bosons at the top of the 
research agenda. Among the experimental realizations are 
systems as diverse as micro-mechanical arrays [1], coupled 
Joscphson Junctions [2], and Bosc-Einstein Condensates 
(BEC) loaded in optical lattices [3], the latter being sug- 
gested as viable candidates for quantum computation. 

Apart from their technological applications, these sys- 
tems offer unique opportunities to experimentally investi- 
gate fundamental problems of quantum mechanics. One of 
the main questions is the manipulation of coherence and 
the stability of complex quantum dynamics under external 
(environmental) perturbations. A well-established mea- 
sure for the latter is the so-called fidelity (also known as 
the Loschmidt Echo). It was first introduced by Peres [4] 
who used fidelity to study quantum-classical correspon- 
dence and identify traces of classical (chaotic or integrable) 
dynamics in quantized systems. Recently, fidelity has 
been used to quantify the corruption of quantum infor- 
mation [5] , in the framework of quantum computation. 

In this Letter, we study the fidelity decay of (ultra-) cold 
atoms loaded in an optical lattice (OL) or coupled micro- 
traps subject to perturbations of the coupling: fco —> ka + 
8k. In the context of OLs this perturbation is achieved by 
adjusting the intensity of the laser beams that create the 
lattice. Such systems arc described by the Bose-Hubbard 



Hamiltonian (BHH), which in second quantization reads 
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where / is the number of sites and denotes summa- 
tion over adjacent sites and U = 4irh 2 a s V ff jm describes 
the interaction between two atoms on a single site (V e g is 
the effective mode volume of each site, m is the atomic 
mass, and a s is the s-wave atomic scattering length). The 
operators hi — b\bi count the number of bosons at site 
i; the annihilation and creation operators hi and h\ obey 
the canonical commutation relations = S^j. In the 

context of Josephson Junction arrays, k is given by the 
Joscphson energy Ej while U accounts for the Coulomb 
interaction of the charged bosons [6]. 

We will focus our presentation on the trimcr / = 3. 
This is the minimum BHH model that contains all the 
generic ingredients [7] (like classical chaotic dynamics [8] ) 
of large BHH lattices and therefore often is used as a proto- 
type model [9] . Our main results arc summarized in fig. 1 
where we plot the fidelity F(t) as function of time for var- 
ious perturbation strengths 5k. The new striking feature 
is the appearance of echoes at multiples of t C cho- By ana- 
lyzing the energy landscape of the perturbation operator, 
we are able to identify i e cho and control the echo efficiency 
by an appropriate choice of the initial preparation. For 
moderate perturbations 8k, an improved RMT modeling 
that incorporates the semiclassical structures of the per- 
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The latter is obtained from Eq. (1) in the limit N 3> 1, 
and reads [13] 




Fig. 1: (Color online) Parametric evolution of the fidelity F(t) 
for different perturbation strengths 8k as a function of time. 
The fidelity exhibits echoes at multiples of t = t cc ho- Here the 
Heisenberg time is tu — 0.55. 



turbation operator can reproduce these echoes, while for 
larger perturbations we rely on semiclassical considera- 
tions. We show that the trajectories leading to fidelity 
echoes becomes more abundant at high energies, in con- 
trast to recent experimental results on echo spectroscopy 
on ultra-cold atoms in atom-optics billiards [10] where it 
was found that the echoes do not survive in the strong per- 
turbation limit. Our analysis indicates that this is due to 
self-trapping phenomena [11] and reflects a generic prop- 
erty of interacting bosons loaded on a lattice. 

The fidelity is defined as the overlap of two states ob- 
tained after time-evolving the same initial preparation 
l^o) with two slightly perturbed Hamiltonians 
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In our analysis, we consider that Hi 2 = Hq =p 5kB. The 
unperturbed Hamiltonian Hq is given by eq. (1) with k = 
fco- In the simulations below we will use fco ~ 15, U ~ 280. 
The number of particles in the simulations is N = 230. 
The coupling operator B is defined as 



B 
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This perturbation is similar to a momentum boost. Such 
a perturbation has recently been investigated [12] in the 
context of the fidelity, where it was found that the fidelity 
freezes at some finite value, as long as the boost is not too 
large. 

Quantum mechanically, we work in the Hq eigenbasis. 
The initial preparation |\I>o) is chosen to be either a sin- 
gle eigenstate of Hq or a Gaussian superposition of them 
centered at energy Eq. A fixed assumption of this work is 
that the perturbation is classically small 6k <C Sk c \, i.e., 
the corresponding classical Hamiltonians Hq, Hi and H2 
are generators of classical dynamics of the same nature. 



H = H/NU 




(4) 



where {Ii, </>j} are the canonical action-angle variables, and 
U = NU. The classical dynamics is governed by the di- 
mensionless ratio A = k/U [8,13]. The other dimensionlcss 
parameter that affects the classical motion is the energy 
E = H. 

An important observable is the generalized force !F(t) = 
-(dH/dX) = J2(tj) \fhXj cxp[i(</?j - tfi)] where t = Ut is 
a rescaled time. We can characterize its fluctuations by 
the autocorrelation function C(i). Its Fourier transform 
C(cl>) is the power spectrum of the fluctuations. Using a 
simple semiclassical recipe, we can connect C(u) with the 
band-profile of the perturbation matrix B [14]: 

a 2 = <|B nm | 2 ) B « (N 2 /U)A ■ C(io nm )/27r (5) 

where we use the basis which is determined by Hq. Above 
(• • • )e indicates an averaging over an energy window SE 
around some energy E, while A ~ U/N [13] is the mean 
level spacing. In cases of chaotic dynamics C(ui) has a 
cutoff frequency lo c = 2t:/t c which is characterized by a 
finite correlation time t c . Equation (5) implies that B is 
a banded matrix with a bandwidth = Huj c . 

Using Linear Response Theory [15, 16] one can obtain 
the following standard expression for the fidelity: 



F(t) ~ 1 - 5k 2 C(t) + 



exp(-Sk 2 C(t)), (6) 



where C(t) = /„' / * C(t' - t")dt'dt". This expression im- 
plies a short time Gaussian decay F(t) = exp[— 8C(0) x 
(8k ■ t) 2 }, which evolves into a long time Gaussian decay 
F(t) = exp[— (a5k) 2 t 2 ] for quantum mechanically small 
perturbations 5k < 5k qm ~ A/cr [16,17,20,21]. The decay 
of F(t) for small perturbations 5k < <5fc qm , is the same 
irrespective of the nature (integrable or chaotic) of the 
underlying classical dynamics. 

If the band-profile of the perturbation operator B were 
flat, we would expect [16, 17, 20, 21] an exponential de- 
cay F(t) ~ exp[— 2Tt] for stronger perturbations <5fc qm < 
5k < Sk plt - U/N. The rate T - (a5k) 2 /A is given 
by the width of the Local Density of States (LDoS) [13]. 
In the perturbative regime (i.e. 5k < 5k prt ), traditional 
RMT considerations [16, 17] can describe the fidelity de- 
cay for chaotic systems [16,17,22]. For even stronger per- 
turbations 5k > Skprt, we enter the semiclassical regime 
[17,20,21]. This identification of the non-perturbative 
regime with the semiclassical limit can be easily under- 
stood if we realize that 6k pi t ~ U/N — > in the limit 
N 1 (while U = const. [13]). Then for any fixed per- 
turbation 5k, eventually 5k 3> <5fc prt . In this regime (and 
provided that the classical dynamics is chaotic) the fidelity 
decay is exponential F(t) ~ exp(— jt), with a rate given 
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Fig. 2: (Color online) The quantum band-profile 
(|B nm | 2 )?7/./V 2 A and classical power spectrum C(uj) (drawn 
with dashed and solid lines respectively) at various eigenvalues 



E { ° ] of the reference Hamiltonian H : a) E% } = 0.22, b) 
p (0) = 0.26, c) E ( n 0) = 0.39. In panels d)-f) we plot the 



E, 

corresponding fidelity F(t) for three representative values 
Sk < 8k a 



< Sk clln (dash-dotted line, deep perturbative regime), 



< Sk < 5k 



pit 



(dashed line, perturbative regime), 
and 8k > 5k pl t (solid line, non-perturbative regime). The 
respective numerical values are for d) and e) 5k = 0.1, 0.5, 2.5, 
and for f) 5k — 0.5, 2.5, 7.5. The initial preparation is an 
eigenstate of Ho. The vertical dashed lines in sub-figures d)-f) 
denote the revival time t ec ho while in a)-c) they denote the 
corresponding frequency uj* — 2n/t ec ho- 



by the classical Lyapunov exponent 7 [23]. Our numeri- 
cal data reported in figs. 2d,e, confirm the validity of the 
above expectations for the BHH model (1) in the chaotic 
regime. 

Less clear is the situation for (predominantly) integrable 
systems for perturbation strengths Sk > Sk prt : depending 
on the initial state, the fidelity decays either faster than 
Gaussian (when the perturbation changes the frequencies 
of the phase-space tori) [15], or in a power law fashion 
(when the primary effect of the perturbation is to change 
the shape of the phase-space tori) [19]. Indeed, our numer- 
ical results for the BHH in the (predominantly) integrable 
regime (high energies) confirm the above expectation (see 
figs. 2f, 3 and related discussion below). 

Next, we focus on the appearance of fidelity echoes. In 
order to gain some insight we turn first to the analysis of 
the perturbation matrix B which generates the dynamics. 
In figs. 2a-c we show the band-profile (|B nm | 2 )£ (dashed 
line) for various energy regimes E. One observes that 
the band-profile is not flat, but exhibits pronounced struc- 
tures within the bandwidth Af,. A striking feature is the 
existence of side-bands whose position w C cho increases as 
we increase the energy E. As a result C(t) (see eq. (6)) 
oscillates, leading to strong fidelity echoes at multiples 
of a characteristic time i e cho = 27r/u; cc i 10 (see figs. 2d-f). 
These echoes arc different from the standard mcsoscopic 



echoes at the Heisenberg time tn = 1/A of quantum sys- 
tems with chaotic classical dynamics [18]. They are instead 
associated with "non-universal" structures that dominate 
the band-profile of the perturbation matrix B and are the 
fingerprints of the lattice confinement. Quantum mechan- 
ically, these are reflected in selection rules that determine 
to which states \m(ko ± Sk)) of Hi 2 an initial state \^o) 
is coupled. 

We consider the case where the initial preparation 
l^o) = \ n (ko)) is an eigenstate of Hq corresponding to 
an eigenvalue E^ ■ Here, the information about the 
coupling of the initial states to the perturbed states 
m(fco ± Sk)), is encoded in the structure of LDoS Pl(E) = 
J2 m \ {m\n)\ 2 S[E - (E m - E n )]. For Sk < Sk prt we have 
that \(m\n)\ 2 « <J„, m + Sk 2 (\B nm \ 2 )/[T 2 + (E n -E m )% 
with a width F « (<rSk) 2 /A < Af, indicating the energy 
regime where most of the probability is contained [13]. For 
short times t < this core will not have time to de- 

phasc/decay. Therefore, the fidelity will show large echo 
recoveries at the time the side-bands have acquired a phase 
of 27r, if t oc ho <C r _1 . The numerical data for Sk < Sk prt 
reported in figs. 2d-f (sec dashed lines) confirm this predic- 
tion. We see that the echo efficiency (i.e. the recovery level 
of F(t)) increases as i cho becomes smaller [24]. Arguing 
along the same lines as above, we are able to explain the 
reduction of fidelity echoes in the non-perturbative regime, 
Sk > Skp T t, observed in figs. 2d-f (see solid lines). Here, the 
LDoS covers the whole bandwidth, spoiling (with respect 
to the perturbative cases) the echo efficiency. 

Similar echoes were experimentally observed in [10] for 
atoms in optical billiards. However, in contrast to our case 
(see black lines in figs. 2d-f), Ref. [10] was reporting a to- 
tal loss of echoes in the semiclassical regime Sk > 5k pi t- 
To understand the origin of the echoes observed for the 
BHH model (1) for Sk > <5fc prt , we will employ semiclas- 
sical considerations. In fig. 3, we compare the quantum 
calculations for Sk > 5k plt with the classical fidelity 

F c i(t) = / dx p- 5k (x,t)p + 8k(x,t) (7) 

where p±su is the classical density function evolved un- 
der H(k ± Sk) and the integration is performed over the 
whole phase space f2 [25]. Although the quantum and 
classical calculations - as far as the echoes are concerned 
- agree quite well, a discrepancy on the fidelity decay be- 
tween F(t) and F c \(t) for small times is also evident. In 
contrast to the chaotic regimes (lower energies), where for 
Sk > 5k pi: t the fidelity decay was exponential with a rate 
given by the classical Lyapunov exponent (see for example 
figs. 2d,e), here F c i(t) ~ t~ a while the quantum fidelity de- 
cays as F(t) ~ i~ 3Q / 2 . The power law decay is a signature 
of classically (predominantly) integrable dynamics while 
the anomalous (faster) quantum decay is a pure quantum 
phenomenon, as was pointed out in [19]. According to 
the prediction of [17], the classical power law exponent 
had to be a = d, where d = 3 is the dimensionality of 



p-3 



Joshua D. Bodyfcltl, Moritz Hillerl,2,3 Tsampikos Kottosl,2 




Fig. 3: (Color online) The quantum (blue line) fidelity decay 
F(t) at E„ — 0.39, plotted together with the correspond- 
ing classical (black line) fidelity decay F c \(t), calculated from 
eq. (7). The perturbation strength is Sk — 7.5 > Sk prt . The 
dashed lines are the best least square fits and are drawn to 
guide the eye. The upper one corresponds to t~ 1,46 while the 
lower one corresponds to t~ ' 99 . 

the system (in fact, in the case of the BHH model (4) 
we have a six-dimensional phase space with two constants 
of motion, namely the total number of particles N and 
the energy E and thus d e g = 2). However, in the high 
energy regime, the dynamics is dictated by self-trapping 
phenomena, leading to localization of particles in one site. 
Therefore the effective dimensionality of the system de- 
scribed by the Hamiltonian (4) is d g = 1. Indeed, the 
best linear fit to the numerical data yields a = 1, thus 
confirming the above argument. 

Let us now return to the analysis of the classical echoes. 
The classical trajectories contributing to the ensemble av- 
erage of the echo, are those that - after evolving forward 
in time with Hi and then backwards in time with Ti.2 - 
return to the vicinity of their initial position. Since Ti.\ 
and Ti.2 arc different in the coupling between nearby wells, 
those classical trajectories that do not jump between wells 
will not feel the difference in the coupling terms. There- 
fore they retrace their forward propagation backwards in 
time, causing the action integral to vanish. These trapped 
trajectories give a perfect contribution to the echo signal. 
Their existence is due to both the discreteness and the 
nonlincarity of the underlying equations of motion. For 
high enough energies the configuration of bosons, local- 
ized in one lattice site persists for long times as can be 
seen from fig. 4. In such cases the bosons are said to be 
"self-trapped". As a result, the echo efficiency is increased 
in agreement with the numerical simulations presented in 
figs. 2d-f. 

We further confirm the importance of the structural av- 
erage band-profile in the appearance of echoes, by eval- 
uating the fidelity decay F(t) using an improved RMT 
(IRMT) model [22]. While traditional RMT models in- 
volve band-profiles with matrix elements given by a Gaus- 





Fig. 4: (Color online) Time-evolution of N = 230 particles 
(bosons) associated with various energy configurations: (top) 
E — 0.06 close to the ground state (equipartition) of the sys- 
tem; (center) E — 0.26 corresponding to a state in the middle 
of the energy spectrum and (bottom) E = 0.39 associated with 
a high energy state. The red line indicates the accessible phase 
space if all particles were on the corresponding site. Note that 
as the energy increases, the likelihood of "trapping" particles, 
i. e. particles that are initially localized in one site (in the spe- 
cific case i = l) and stay there for long times, increases. 

sian distribution with a constant variance a 2 within the 
bandwidth [16, 18], the IRMT model incorporates the 
band-structure associated with the dynamical system, 
through the semiclassical relation (5). In fig. 5 we present 
the results of these simulations. An excellent agreement 
with the fidelity calculations of the BHH (1) for the pertur- 
bativc regime Sk < Sk pi t is evident. At the same time the 
IRMT modeling cannot describe the quantum results in 
the non-perturbative regime. Here, the echoes are related 
to subtler correlations of dynamical nature (self-trapping), 
which goes beyond the autocorrelation function C (r) that 
determines the band-profile. 

Finally, in fig. 5 we report our numerical calculations 
for the case where \^o) is a superposition of eigenstates 
of Hq with a Gaussian weight around some energy Eq. 
We found that as the width a 2 G of the Gaussian prepara- 
tion increases, such that the participating eigenstates of 
Hq span an energy window which is wider than SE over 
which the band-profile of B is constant (see for example 
Figs. 2a-c), the echo-efficiency is suppressed. This is due 
to destructive interferences resulting from states located 
at various parts of the spectrum where the band-structure 
(thus the w cc ho) is different. The destruction of echo effi- 
ciency becomes more pronounced in the non-perturbative 
regime where no sign of echo is observed (see fig. 5). 

Although the above analysis is focused on the / = 3 
case, the reported results are expected to be applicable 
for larger lattices / > 3 as well. Furthermore, due to the 
connection between fidelity decay and loss of coherence 
[27], our findings can be used to engineer coherence echoes 
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Fig. 5: (Color online) The fidelity F(t) at E^ = 0.26 for two 
perturbation strengths (upper curves 5k = 0.5 < Sk prt , lower 
curves 8k — 2.5 > <5fc prt ) and various initial conditions: black 
solid lines correspond to 8 states, dashed lines to Gaussian 
states with variance a 2 = 10, and dash-dotted lines to 5-states 
evolved under the IRMT model. We have found that within 
the perturbative regime, the IRMT calculations coincide with 
the Linear Response Theory expression from eq. (6). 



of a central system coupled to a quantum bath consisting 
of cold atoms in an optical lattice [10,26]. 
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